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!PREPROCESADO, ENTRADA DE DATOS 

prep7 

window. l.full  ! define el numero v tamaño de las 
! ventanas gráficas 

view,1,0,0.9 

vup,l,y ' especifica la orientación de los ejes globales 
! eje Y vertical 

title, MURO PANTALLA 


'TIPO DE ELEMENTO 


et,1,42 ! elemento tipo plane42 
! elemento plano con 4 nodos 
kevopt,1,3.3 ! sección tipo 1. opción 3, valor 3 
! keyopt3=3 tensiones planas 
r,1,250 ' espesor del elemento [mm]. 


! valor requerido si kevopt3=3 


¡MATERIAL hormigón 


mp,ex 1,20 ! modulo de elasticidad [kKN:mm2] material 1 

mp,nuxy,1,.1  !coef. de poisson material 1 

mp,dens, 1,0 |! densidad: no considera el peso propio 
¡NODOS 


n,1,0,36000.0  ! coordenadas del nodo 1 en (0,36000.0) 
ngen,9,1,1,1,1,750,0,0 ! generación de nodos 
ngen,37,9,1,9.1,0.-1000,0 


"ELEMENTOS 
type, l l asigna el tipo | a los siguientes elementos 
real, l 
e.10,11,2,1 ! define el elemento entre los nodos 
110,11,2,y 1 
egen,8,1,1,1,1  ! generación de elementos 
egen,3,9,1,2.1 
egen,3,9,7,8.1 
egen,12,27.1.16.1 


¡CONDICIONES DE APOYO 
nsel,s.loc,v.0 


d,allux.0 
d.all.uv,0 
nsel, all 
"CARGAS 
nsel,s.loc.x.0 ' selecciona el subconjunto (s) de nodos 
' con coordenada local x= 0 
Lall,Éx,7.5 ' aplica cargas puntuales en (all) todos los nodos 


! seleccionados en el sentido del eje X 


ESTRUCTURAS II] 


! de valor 7.5 [KN] 
nsel.all 
nsel.s,loc,x,6000 
fall,£x.7.5 
nsel,all 


¡SALVAR Y SALIR 
save 
finish ! salir del preprocesador 


TIPO DE ELEMENTO 
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"PROCESADOR DE SOLUCIÓN 
solution 
antvpc. static | indica que el tipo de análisis es estático 
solve ' ordena resolver el sistema 
finish ! salir del procesador 


'POSTPROCESADOR 

'SALIDA E INTERPRETACIÓN DE RESULTADOS 
post] 
etable.sx.s,x ! prepara la tabla denominada sx. 
! con los valores de los esfuerzos s componente x 
etable.sv.s.v ! ídem. componente y 


etable.sxv.s.xv  ! ídem componente xv 


set.!.l ! lee datos de resultados del paso de carga 1 
'DIBUJAR 
pldisp,2 ! muestra la estructura deformada 


! con Ja estructura original de base (2) 

! gráfica de los datos de la tabla sx. 

! (avg): tomando en el interior del elemento un 

! valor como promedio lineal 

plctab.sv.avg  ! idem tabla ss 

pletab.sxv.avg  ! idem tabla sx» 

plvect.pdir ! gráfica de resultados del elemento como vector 
! (pdir): tensiones en direcciones principales 

! gráfica de los desplazamientos de los nodos 

! como vector 


pletab.sxavg 


plvect.disp 


'IMPRIMIR 
pretab.sx.Sy.SXY ! imprime los datos de las tablas sx. sy. 
'sxv [N mm] 
' imprime las reacciones en los nudos, (f) fuerzas 
! imprime los esfuerzos de los elementos 
! (f): fuerzas en el sentido de los ejes X, Y. Z 


finish ! salir del postprocesador 
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CORTADURA MÁXIMA 7 
Tmox= 1. 1N/mm* 


se produce en los dinteles de los niveles 1 y 2 


zona 1 
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Agrupamos las fuerzas nodales de las secciones 
en las zonas 1 y 2 en una resultante 


ZONA 1 


Con estos esfuerzos podemos tener una idea en 
relacion a los modelos a, b y c. 
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CONJETURAS: 
TRES MODELOS DE COMPRESIÓN 


a) sin huecos 


MOMENTOS EN LOS SOPORTES 
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Ma=9720 kN-m 
20kKN mm 


o TENSIONES DEFORMACIÓN 
= 15 m:3Em” - 9720kN.m lg 2 y 29720 KN.m. A 36m 


e 
: 6 m 250mm 
12 


V max 


9720kN.m 


TO IET DÓ 3N má 
6000” *0256mnrim b) Como 2 ménsulas la deformación será 16 veces la anterior 


3=736mm 


| c) Estimamos la deformación por cortante 
>. LL AL TL? 
lo i — L 1El 
1 : | ; > 
le | o 
E pación 


Mb=1/2 Ma AAA =014mm/m 
12-20 kN /mm' S-1250. pop! 


rap, 
Pb=- 3 Wa =5¿HWa y 
0,=8- 0, =8:6,5=52N/mm” as. incremento de deformación: 
AS =0,14mm.m+*36m.2= 2,6 mm 


| a, | 
CORTADURA 
720kN.m+ 13 
Ñ xd D) y ION me 36m/2 _ 829 3 
o ; 
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ESTRUCTURAS III ANALISIS MANUAL 


8. MURO PANTALLA 


4.42 PLANE42 2-D Structural Solid 


Table 4.42-1 PLANE42 


Element Name 
Nodes 
Degrees ofFreedom 


Real Constants 


Material Properties 


K,L 


Surface Loads 


Element Coordinate 


System (shown for 
(or Axial) KEYOPT(1)=1) J Body Loads 
X 


| I 
(Trianguler Option - not recommended) 
J 
X (or Radial) (1) 


Figure 4.42-1 PLANE42 2-D Structural Solid Special Features 
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Input Summary 


PLANE42 
LLL 
UX, UY 


None, iFKEYOPT(2) = 0, 1, 2 
Thickness, if KEYOPT(S) =3 


EX, EY, EZ, (PRXY, PRYZ, PRXZ or NUXY, NUYZ, 
NUXZ), ALPX, ALPY ALPZ, DENS, GXY, DAMP 


Pressures: 


face 1 (J-1), face 2 (K-J), face 3 (L-K), face 4 
(I-L) 


Temperatures: 

11), TJ), TK), TL) 
Fluences: 

FL(1), PLC), FLECK), FLEL) 


Plasticity, Creep, Swelling, Stress stiffening, Large 
deflection, Large strain, Birth and death, Adaptive 
descent 
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!¡PREPROCESADO, ENTRADA DE DATOS 
/prep7 
/view,1,2,3,5 
/vup,1,z 
/title, LOSA ANFITEATRO 

'TIPO DE ELEMENTO 
et,1,63 ! elemento tipo shell63. Comportamiento de membrana 

! y de flexión con 4 nodos 

1,1,450 ! espesor [mm] 

! MATERIAL 

! asimilación a hormigon armado de cuantia alta 

mp,ex,1,21 ! [kN/mm?] 
mp,nuxy,1,.1 ] 
mp,dens,1,8e-9 ! [kN/mm*] 

! GEOMETRIA 
csys,1 ! activa el sistema de coordenadas cilindricas (1) 

!puntos clave 
k,1,4000,45 ! punto clave numero 1 de coordenadas cilindricas 

! (R=4000, Angulo=45) [mm] 

k,2,15000,45 
k,3,15000,-90 
k,4,4000,-90 

lineas 
1,1,2,6 ! define una línea entre los puntos clave 1 y 2 
idividiéndola en 6 partes; esta será la línea 1 
1:23:12 
1,3,4,6 
1,4,1,12 
Iplot ! dibuja las lineas definidas 
al,1,2,3,4  ! define un area delimitada por las lineas 1, 2, 3 y 4 

! sera el area 1 

! GENERACION DE NODOS Y ELEMENTOS 
type, 1 
real, 1 
amesh,1 ! genera elementos superficiales con sus nodos 

! en el area 1, 
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O TT 


! asignando el elemento seleccionado 
/pnum,elem,]1 ! activa la numeracion de elementos 
/pnum,node,1  ! activca la numeracion de nodos 
eplot ! dibuja los elementos 


! Condiciones de apoyo 
/pbc,all,1 ! muestra los simbolos de coacciones 
nsel,s,loc,x,15000 
d,all,uz,0 í apoyo continuo en el contorno 
nsel,all 
d,1,uz,0 
nsel,s,loc,y,-90 
d,all,ux,0 | se impide movimientos horizontales (x) 
d,all,roty,0 ! y giros por simetria 
ldsym,symm,x !también se puede decir esta orden en vez de 
!las dos anteriores (plano de simetria perpendicular al eje x) 
d,8,uy,0 limpide mov. Horiz. Eje y 
nsel,all . 
¡CARGAS 
! losa aligerada 8e-9kN/mm' * 450mm= 3.6e-6kN/mm? 
! sobrecarga 4e-6 kN/mm? 
t 
! 7.6e-6kN/mm? 
! o bien 2.11*peso propio 
acel,0,0,2.11 ! asigna una aceleracion 2.11 en el eje Z 
Isfe,all,1,pres,,7.6e-6 lotra forma de introducir la carga 
! carga lineal 15.2e-3kN/mm 
nsel,s,loc,x,4000  ! selecciona los nodos de radio = 4000mm 
! los del perimetro interior 
l aplica una carga puntual en el sentido del 
! eje Z en cada uno de ellos de -11.94 [kN] 


f,all,fz,-11.94 


nsel,all 
ISALVAR Y SALIR 
save 


finish ! salir del preprocesador 
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¡PROCESADOR DE SOLUCION MOMENTOS RADIALES M, DEFORMACIONES Y REACCIONES 
/solution 

antype, static ! indica que el tipo de análisis es estático 

solve ! ordena resolver el sistema 

finish ! salir del procesador 


'POSTPROCESADOR 
ISALIDA E INTERPRETACION DE RESULTADOS 
/postl 
! Tablas de valores 
etable,sx,s,X ! prepara la tabla denominada sx, 
! con los valores de los esfuerzos s componente x 
! en coordenadas globales 
etable,sy,S,y 
etable,sxy,S,Xy 
set,1,1 ! lee datos de resultados del paso de carga 1 
smult,mx,sx,,33750 
! calcula los momentos a partir de las tensiones 
IM=S*W 
! en coordenadas globales 
smult,my,Sy,,33750 
smult,mxy,sxy,,33750 


etable,m1,smisc,4 NL 
! momento en coordenadas locales del elemento DIRECCIONES PRINCIPALES 


etable,m2,smisc,5 zona 1: moemntos de igual signo 
etable,m12,smisc,6 o zona 2 momentos de signo distinto 
! Dibujos 
pldisp,2 ! muestra la estructura deformada 
! y el contorno de la original (2) 
/pbc,al1,0 ! no muestra los simbolos de coacciones 
pletab,mx,avg  ! gráfica de los datos de la tabla mx, 


! (avg): tomando en el interior del elemento MOMENTOS CRUZADOS My 


! un promedio lineal Me 
pletab,my,avg Msg, ma 239 KN.m/m 


pletab,mxy,avg 

pletab,m1,avg 

pletab,m2,avg 

pletab,m12,avg 

plvect,pdir ! gráfica de resultados del elemento como vector 

! (pdir): tensiones en direcciones principales 

limprimir 

pretab,mx,my,mxy 

pretab,m1,m2,m12 

prrsol,f ! imprime las reacciones en los nudos, (£) 

prnsol,u ! imprime desplazamientos (u) en los nodos 

finish 
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9. LOSA ANFITEATRO 9.2 


ro AR APOYADA EN CONTORNO 3) LOSA CON AGUJERO Y CARGA LINEAL 5) ASIMILACION A ESTRUCTURA DE VIGAS 
Solución elástica (v=0) a) vigas radiales con anillo 
19. ns T,= ro/r:q Las vigas se apoyan ficticiamente en 
=p ESRP+SR) M,= -1yq/2:(Ln r+1,5)+C,/2-C,/r2 apuyos Pvc ás 
M.= -19'q/2-(Ln r-0,5)+C,/2+C,/r? distancia. (Si tomamos máxima 15m) 


2 
5 MISS 130m/m 


:q/16:(R2-r2) “Los pilares soportan más de la mitad de 
a= q/16:(3R2-r2) 8 la carga cuando se aproximan al eje de 
: simetría. Y el central menos de Y de 
¡ carga (p. ej 0,4) 
2 KA 00 fan, ' 


m/m 


CARGA TOTAL 4029 kN 
Y, CARGA TOTAL 2013 kN PILAR CENTRAL 805 kN 


2) LOSA CON AGUJERO 4) LOSA CON APOYO INTERIOR 


b) vigas en luces cortas 
Solución elástica 


Llegamos de esta manera de nuevo a 
T,= q:1/2-q:1? T,= q:1/2-q:1,219:N/r luces de flexión de 15m. 


¡= -3q:12/16+q:1p2/4-(Ln r+0,5)+C,/2-Cy/r2  M,= -3-q:r2/16+q'rp2/4:(Ln r+0,5) +r9N/2-(Ln r+1,5)+C,/2-C,/r2 
a= -q:12/16+q:rp2/4-(Ln r-0,5)+C,/2+C,/r2 M.= -q:1/16+q:r,2/4:(Ln r-0,5)+r,N/2-(Ln r-0,5)+C,/2+C,/r2 


259 K£ mfun 
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9. LOSA ANFITEATRO 9.3 


¡PREPROCESADO, ENTRADA DE DATOS 
/prep7 

/view,1,2,3,5 

/vup,1,z 

/title, LOSA CON AGUJERO Y DOBLE APOYO 


¡TIPO DE ELEMENTO 
et 1,63  ! elemento tipo shell63 
! comportamiento de membrana y de flexión 
! con 4 nodos 
r,1,450  ! espesor [mm] 


! MATERIAL 
! asimilacion a hormigon armado de cuantia alta 

mp,ex,1,21 ! [KN/mm2] 
mp,nuxy,1,.1 
mp,dens,1,8e-9 ! [KN/mm3] 

! GEOMETRIA 
csys,1 l activa el sistema de coordenadas cilindricas 
(1) 

Ipuntos clave 
k,1,4000,90 ! punto clave numero 1 de coordenadas 
cilindricas (4000,45) [mm] 
k,2,15000,90 
k,3,15000,0 
k,4,4000,0 

lineas 
1,1,2,20  ! define la linea 1 entre los puntos clave 2 y 6 
1,2,3,9 
1,3,4,20 
1,4,1,9 
Iplot ! dibuja las lineas definidas 
al,1,2,3,4 ! define un area delimitada por las lineas 1, 2, 3 y 4 

! sera el area 1 


! GENERACION DE NODOS Y ELEMENTOS 

type, 1 

real,1 

amesh,l ! genera elementos superficiales con sus nodos 
! enel area 1, 
! asignando el elemento seleccionado 

/pnum,elem,1 ! activa la numeracion de elementos 

/pnum,node,1 l activca la numeracion de nodos 

eplot ! dibuja los elementos 


! Condiciones de apoyo 
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/pbc,all, 1 l muestra los simbolos de coacciones 
nsel,s,loc,x,15000 

d,all,uz,O ! apoyo continuo en el contorno 
nsel,all 

nsel,s,loc,x,4000 

d,all,uz,0 

nsel,all 

nsel,s,loc,y,90 

dsym,symm,x 

nsel,all 

nsel,s,loc,y,0 

dsym,Symm, y 

nsel,all 


¡CARGAS 
i losa aligerada 8e-9kN/mm3 * 450mm =3.6e-6kN/mm2 
! sobrecarga 
4e-6 kN/mm2 


7.6e-6kN/mm2 

! 

o bien 2.11*peso propio 
acel,0,0,2.11 ! asigna una aceleracion 2.11 en el eje Z 
nsel,all 


ISALVAR Y SALIR 
save 
finish ! salir del preprocesador 


¡PROCESADOR DE SOLUCION 
/solution 
antype, static l indica que el tipo de análisis es estático 
solve ! ordena resolver el sistema 
finish ! salir del procesador 


¡POSTPROCESADOR 

ISALIDA E INTERPRETACION DE RESULTADOS 
/postl 

! Tablas de valores 
etable,m1,smisc,4 


! momento en coordenadas locales del elemento 


etable,m2,smisc,5 
etable,m12,smisc,6 
! Dibujos 


pidisp,2 ! muestra la estructura deformada 

! y el contorno de la original (2) 
/pbc,al1,0 ! no muestra los simbolos de coacciones 
pletab,m1,avg 
pletab,m2,avg 
pletab,m12,avg 
plvect,pdir ! gráfica de resultados del elemento como vector 
! (pdir): tensiones en direcciones principales 

limprimir 

pretab,m1,m2,m12 
prrsol,f ! imprime las reacciones en los nudos, (f) 
prnsol,u ! imprime desplazamientos (u) en los nodos 


¡OTRAS ORDENES PARA SACAR GRAFICAS 
/erase 

/yrange,-55,150 
Ipath,1,2 

pdef,cero,s,z 
plpath,cero 

/moerase 
etable,m1,smisc,4 
Ipath, 1,2 
pdef,mrad,etab,m1,avg 
plpath,mrad 


/erase 

/yrange,0,160 
pdef,cero,s,z 
plpath,cero 

/noerase 
etable,m2,smisc,5 
Ipath, 1,2 
pdef,mcirc,etab,m2,avg 
plpath,mcirc 


455 ESCUELA TÉCNICA SUPERIOR DE ARQUITECTURA 
U.P.M. MADRID 


ESTRUCTURAS II EN DE LOS ELEMENTOS FINITOS 


9. LOSA ANFITEATRO 9.4 


IGLESIA DE LA ATLÁNTIDA. 


ELADIO DIESTE 


URUGUAY. 1960 


CORTE TRANSVERSAL 


WAVE 
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DECIR 


“Puepo ÉSTE CEL PRoyecTo) RESPONDIERA A 


QuE PRocorRÉ QuE 


UN ESTILO SERIO, A LA VER SEVERO Y AMABLE DE PIEPAD cow UNA 


GRACO CONFIANZA En UL ESPIRITU CRISTIANO DE LOS MUMILOES QuE 


HAN DE USARLA .., Que LA IGLÉSIA (omo ARQUITECTORA ND FUEST 


UL ORSTA COLO PARA UA PIEDAD VERDADERA 5D dv MANIFESTA Cl0J 


PRIMERA LoS MEDIOS y los MATERIALES USADOS EL LA CONSTROCC ¡Sn 


de Cm) 


pl, . 
TAMBIEN QUIEREN SER. EXPRESJIOS: HUMILDES COMO Los PieLes 


PARA QUIENES LA IGLEJIA Sie COMSTRUYy” 
2 


PERO FUERO TRATA-DOS 


low El OÉSVE Lo  Gue ASPIRA AÁSER El HoMENAJE QUÉ ESTOS 


HUMILDES s€ MERECEN ds 
¡A Ass Dieske 
METODO DE LOS ELEMENTOS FINITOS 
ESTRUCTURAS III pri ] 
q 10. IGLESIA DE LA ATLANTIDA 


10.1 


! PREPROCESADOR: Entrada de datos 
/prep7 
/window, 1,full 
/view,1,20,20,20 
/vup,l.z 
/title, LA ATLÁNTIDA 


! TIPO DE ELEMENTO 
et,1,63 ! elemento tipo shel163 
! plano de 4 nodos 
r,1,0.1 l espesor cubierta [m] 
r,2,0.3 ! espesor muro [m] 


! MATERIAL 1: ladrillo 
mp,ex, 1.686 ! módulo de elasticidad 
[kN/m2] 
mp.nuxy,1,e5  ! coef. de Poisson 
mp,dens,1,16  ! densidad [kN/m3] 


! GEOMETRIA 

! definición del muro 

! línea inferior 
k,1,8,-3,0 ! puntos clave 
k,2,8,0,0 


1,1,2 ! L1 recta entre los puntos 1 y 2 


! línea superior 
k.3,6.5,-3,7 
k,4,6.8,-2.25,7 
k,5,8,-1.5,7 
k,6.9.125,-0.75,7 
k,7,9.5,0,7 
bspline,3,4,5,6,7  'L2 curva de segundo 
grado 
l entre los puntos 3,4,5,6,7 
! líneas interiores 
1,13 1L3 
127 1 L4 
l áreas muro 
al,1,4,2,3 1 Al definición de área 
entre 
las líneas 1, 4, 2, 3 
aplot 


! Definición de la cubierta 

! Tínea central 
k,8.0,-3,7.07 Ipuntos clave 
k,9,0,-2.25,7.25 
k,10,0,-1.5,7.787 
k.11.0,-0.75,8.32 
k,12,0,0,8.5 
bspline,8,9,10,11,12 

|! L5 curva de segundo grado 

! parábolas transversales 
k,13,4,-3,7.0 
k,14,4,0,8.236 
bspline,3,13,8 !L6 
bspline,7,14,12 1L7 

l áreas cubierta 
al,2,7,5.6 1 A2 
aplot 


! Definición del zuncho 
! línea exterior 
k,15,7.9.-3,7 
k,16,8.2.-2.25.7 
k,17,8.95,-1.5.7 
k,18,9.6,-0.75,7 
k,19,10,0,7 
bspline,15,16,17,18,19 
! líneas muro-zuncho 
1,3,15 
1.7,19 
l áreas zuncho 
al,8,10,2,9 
aplot 


! Simetrías 

arsym,x,1,3,1,.0 Isimetría respecto al eje x 
! de las areas de la 1 a la 3 tomadas 
! de 1 en 1, generando áreas 
! adicionales 

arsym.,y.1,6,1.,0 Isimetría respecto al eje y 
! de las áreas de la 1 a la 6 tomadas 
! de 1 en 1, generando áreas 
! adicionales 

nummrg,all ! elimina puntos y rectas muy 
! próximos o dobles 
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! NODOS Y -ELEMENTOS 
! muro 
type. 1 
real,2 
! eshape,0,2 ! rellenar el área permitiendo 
! mezclar elementos cuadrados 
y 
! triangulares sin crear nodos 
! intermedios 
l esize,,5 ' número de divisiones de las 
! líneas de borde 
amesh,1,10,3  ! rellena de elementos del 
l área l ala 10 en pasos de 3 


! bóveda 

type, 

real, 

! eshape,2,2 ! sólo permite elementos 
' cuadrados sin nodos 
! intermedios 

l esize,,7 

amesh,2,12,3 


! TIPO DE ELEMENTO 
et,2,63 lelemento tipo shell63 
r,3,0.2 lespesor zuncho [m] 
! MATERIAL 2: hormigón 
mp,ex,2,21e6  !módulo de elasticidad [kN/m2] 
mp,nuxy,2,e5  !coef. de Poisson 
mp,dens,2,25  !densidad [kN/m3 ] 


! zuncho 
type,2 
real,3 
leshape,2,2 
lesize,.3 
amesh,3,12,3 


! CONDICIONES DE APOYO 
! empotramiento en la base 
nsel,z,0 Iselecciona los nodos de z=0 
d,all,uz,0 Icoarta desplazamiento vertical 
d,all,uy,0  !coarta desplazamiento en la 
¡direccion y 
d,all,ux,0 Icoarta desplazamiento en la 
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¡direccion x 
d,all,rotx,0 Icoarta giro respecto eje x 
d,all,roty,0 Icoarta giro respecto eje y 


Sobrecarga a un lado de la estructura 
Carga superficial presión en el eje Z 


Viento: Se supone que la acción horizontal la 
absorbe el zuncho al ser más rígido. 


MALLA DS ELEMENTOS 


+ En AUSy3 5.3 educac. 


45 necesafio 


Bomitar al mm o 


de. edementos 
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ESTRUCTURAS II 


ANSYS 


d,all,rotz,0 Icoarta giro respecto eje z 
nsel,all Iselecciona todos los nodos 


! simetría en planos Y 
nsel, y,-3 Iselecciona los nodos de y=-3 
dsym.symm.y  !coacciones equivalentes a que 
la ¡estructura es simétrica respecto el 

¡eje y 

nsel,y,3 
dsym,symm, y 
nsel,all 


! CARGAS 
! Peso propio 
acel.0,0,1 | carga de peso propio 
! Sobrecarga nieve a un lado 
! sfa,2,.pres,-0.5 
! sfa,8,.pres,-0.5 
! Viento 
! Cargas en nodos del zuncho 
1 f,915,fx.-2.1 
1 £,917,£x,-2.1 
1 f,918,fx,-2.1 
1 f,919.fx.-2.1 
1 f,920,fx,-2.1 
! 


1 £.989,fx,-2.1 
1 £,985,fx,-2.1 


save 
finish 


'PROCESADOR DE SOLUCIÓN 
/solution 
lumpm.on 
antype.static 
solve 

finish 


lanálisis estático 
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'POSTPROCESADOR A 
ISALIDA E INTERPRETACIÓN DE RESULTADOS 


? 
/post!l ; A x 
prrsol.f ! imprime reacciones en nudos - AR 
|! comprobación de equilibrio global =p 2 e a - En 
pldisp.2 ! estructura deformada y contorno de la original (2) +++ BAJE HURO 
prnsol,u ! imprime desplazamientos (u) en nudos 34, 2 DISTRIBUCIÓN De ¡A 
! comprobación de deformaciones ) Peso 
! Tensiones 2), lego los dei ' N 
shell,top ! localización en que se dan los resultados de las uan a A AR A EMMA 16? 
! tensiones en el elemento laminar a e yz0'y kN 
! (top): cara superior Res muro =2x (10 o x 031 4x4)* 
! (mid): medio Pos rado: De As KM)» xo/10,77%7) 5 32,50 
! (bot): cara inferior Tr le +0!5 ku] 
etable.SXt,s,x OA 


etable,S Yt.s.y 
etable,SZT,s.z 


shel!l,bot 

etable,SXb,s,x 
etable,S Yb,s,y 
etable,SZb,s,z 


pletab,SXt,avg 
pletab,SXb,avg 


pletab,SYt.avg 
pletab,SYb,avg 


pletab,SZt.avg 
pletab,SZb,avg 


! Momentos en coordenadas locales del elemento 
etable,Mx,smisc,4 
etable,My,smisc,5 
etable,Mxy,smisc,6 
pletab,Mx,avg 
pletab,My.avg 
pletab,Mxy,avg 
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ER ASS 
plvect.s..,,Vect,,on Tihxls-+ lerase 
113 AS J /yrange.-850.850 
! GRÁFICAS ELIT ATA shell.top 
SE P or ttable,SXt,s,x 

! Momentos y tensiones en cortes transversales Edda et 4 ño shell,mid 
Jerase EN sel E $ y *  etable.SXm.s.x 
/yrange,-1.5,1.5 e shell,bot 
etable,Mx,smisc,4 E A Rs d etable.SXb.s.x 
Ipath,487,488,489,490,491,492.493,472,728,727,7 e + Ipath,487,488,489,490,491,492.493,472,728,727.726.72 
26.725,724.723,708 e 5.724.723.708 
pdef.Mx1,etab,Mx.avg - pdef.SXtop.etab,SXt.avg 
plpath,Mx1 plpath.SXtop 
/noerase Ñ /noerase 
etable.Mx,smisc,4 ; y pdef.SXmid.etab,SXm,avg 
Ipath.498,558,583,572,542.481,767,743,739,751,7 LA TEMSIONES PRINCIPALES plpath.SXmid 


80,713 


LA ATLÁFTIDA 
CORTE TRANSVERSAL CENTRAL 


pdef,SXbot,etab,SXb,avg 


de x,av UN (7 ER Po e E GE] Í 

Er prisa :AVg a Y5 KN/po? ANSYS| — plpath.SXbot 
/noerase A [850 ., /erase 
etable,Mx,smisc,4 2 pe Mx bogoe “o ll Sa her ¡yrange,-850,850 
Ipath,596,671,683,684,638,665,615,861,838,889,8 EN o de a mE ES 2 sxhio shell.top 
88.887,816 O! sx pes etable,SXt.s.x 

qn 5 NA CEnRO ' á 
pdef,Mx3,etab.Mx.avg $2 á shell.mid 
pipath,Mx3 A, 0 etable,SXm.s,x 

di shell,bot 
! Momentos y tensiones en cortes transversales = etable.SXb.s.x 
lerase . Ipath.498.558.583,572,542.481.767.743.739.751.780.71 
/yrange.-1.5,1.5 US TA a Li an 3 
etable.My,smisc,5 EEN DIST | pdef.SXtop.etab.SXt.avg 
Ipath,487,488,489.490,491,492.493,472,728,727,7 > TESTE TRANENERAL INTERMEDIO plpath.SXtop 
26.725,724,723.708 Ds BA AOS cid /noerase 
pdef,My1,etab,My,avg ds ARSS] | KN/m2 ANSYS| — pdef.SXmid.etab,SXm,avg 
plpath,My1l 5 plpath.SXmid 
/noerase A pdef.SXbot,etab.SXb,avg 
etable,My,smisc,5 dá plpath,SXbot 
Ipath,498,558,583,572,542,481,767,743,739,751,7 4 
80,713 | lerase pdef.SXtop.etab.SXt.avg 
pdef.My2,etab,My.avg 5 /yrange,-850,850 plpath,SXtop 
plpath.My2 shell.top Inoerase 
/noerase m. ¿etable,SXt.s.x pdef,S Xmid,etab.SXm.avg 
etable,My,smisc,5 > shell,mid plpath,SXmid 
Ipath,596,671,683,684,638.665.615.861,838,889.8 le PPP. etable,SXm,s.x pdef,SXbot.etab.SXb.avg 
88,887,816 “A 0> A A A A E di pa 0 shell.bot plpath.SXbot 
pdef,My3,etab,My.avg acia el LA ATLANTA eS etable,SXb,s,x 
plpath.My3 My cocos: desmsuersles CORTE TRANINERSAL Ipath,596.671.683.684.638.665.615,861,838,889.888,88 


CERCANO A BORDE 


7,816 
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10.5 


¡Momentos y tensiones en cortes longitudinales 
/erase 

/yrange.-2,2 

etable,Mx.smisc,4 
Ipath.5,473,475,477,479,481.483,485.472,621,619, 
617,615,613,611,609,136 

pdef.Mx4.etab,Mx.avg 

plpath,Mx4 

/noerase 

etable,Mx,smisc,4 
Ipath,709,711,713,715,717,719,721,708,827,825,8 
23,821,819,817,815 

pdef,Mx5,etab,Mx,avg 

plpath.Mx5 


/erase 

/yrange,-2,2 

etable,My,smisc,5 
Ipath,5,473,475,477,479,481,483.485,472,621,619, 
617,615,613,611,609,136 

pdef,My4,etab,My,avg 

plpath,My4 

/noerase 

etable,My,smisc,5 
Ipath,709,711,713,715,717,719,721,708,827,825,8 
23,821,819,817,815 

pdef.MyS,etab,My,avg 

plpath,My5 


shell,mid 

lerase 

/yrange,-850,850 

etable,SX,s,x 
Ipath,5,473,475,477,479.481,483.485,472,621,619, 
617,615,613,611,609.136 

pdef,SX4,etab,SX.avg 

plpath,SX4 

/noerase 

etable.SX,s,x 
Ipath,709,711,713,715,717,719,721,708,827,825,8 
23,821,819,817,815 

pdef,SX5,etab,SX,avg 

plpath,SX5 


2 a 


LA ATLÁNTIDA 


Me CORTE LONGITUDINAL 


i 


knNi.m 


NA 


ZO 


CORTE  LONSITUDINAL 


A NODOS 2ZUYNCMO BÓVEDA 


/erase 

/yrange.-850.850 

etable,SY s.y 
Ipath,5,473,475,477,479.481.483.485,472.621.619.617.6 
15.613.611,609,136 

pdef.SY4.etab,S Y avg 

plpath.SY4 

/noerase 

etable,S Ys, y 
Ipath,709,711,713,715,717,719,721,708.827,825,823,82 
1,819,817,815 

pdef,SY5,etab,SY avg 

plpath,SY5 


!Empuje en zuncho 

erase 

/yrange,-750,200 
Ipath,20,22,24,26,28,30,32,34,5,276,274,272,270,268,26 
6,264,252 

pdef.empuje,s,x 

plpath,empuje 


lErrores 
/erase 
etable,errorl.serr  ! porcentaje de error debido a la 
|! geometría 
etable,error2.sdsg  ! variación máxima de las tensiones 
pletab.errorl avg 
pletab.error2.avg 
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ESTRUCTURAS III 


'DEFINICIÓN APROXIMADA DE LA CÚPULA 


SAN PEDRO 


Figure 4.65-2 


DE SAN PEDRO ¿ 


34, 


HOW (deso fol lnierua) 


¡PREPROCESADO, ENTRADA DE DATOS 


/title, SAN PEDRO 


/view,1,0,1,1 
/vup,1,z 
/prep7 


'TIPO DE ELEMENTO 


et,1 


65 ! SOLID6S, 


l elemento con capacidad de fisuras en tracción 


IMATERIAL ladrillo 


Lo 


mp.ex,1,5e6 
"mp,dens,1,20  ! densidad 20 kN/m3 


' modulo de elasticidad [kN/m2] 


mp,nuxy,1,.15 !nu=0,15 


tb,concr,1 


' tipo hormigon SOLID6S 


tbdata,1,0,1,355,30000 


SOLID65 Stress Output 
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! define los valores de la tabla de datos 

! tabla para fisuracion: material 1, 

! tranferencia de cortante O (fisura abierta), 
1 1 (fisura cerrada), 

! tensión de rotura en tracion 355 [kN/m2], 


U.P.M. MADRID 


l activa tabla de datos del material 1, 


ESTRUCTURAS HI 


! en compresion 30000 [kN/m2] 


¡GEOMETRIA 
¡NODOS 

local,11,2,-3,0,0 ! coordenadas esféricas (2) para 

lentrada de sección 

n,1,19,0,0 

n,21,19,0,70 

£111,1,21 

n,22,19.75,0,0 

n,42,19.5,0,70 

f111,22,42 

n,43,20.5,0,0 

n,63,20,0,70 

£111,43,63 

n,64,21.25,0,0 

n,84,20.5,0,70 

f111,64,84 

n,85,22,0,0 

n,105,21,0,70 

111,85,105 

local,12,1,0,0,0 ! define coordenadas cilindricas (1) 
ipara generar por rotación un octavo 
! de cúpula 

ngen,3,105,1,105,1,0,1,0 ! generando 3 bandas 
! (1 grados de separación entre ellas) 


¡ELEMENTOS 
lesys,0 ! selecciona la orientación de ejes en los 
elementos 
! de cara a la salida de resultados 
ipor defecto (0):paralelas a ejes de 
elementos 
esys,11 
e,1,22,127,106,2,23,128,107 
egen,4,21,1,1 ! 4 elementos en el espesor 
egen,20,1,1,4,1 ! 20 elementos en toda la altura 
egen,2,105,1,80,1 12 elementos 


¡CONDICIONES DE APOYO Y CARGA 

csys,12 l activa el sistema de coordenadas 
cilíndricas previamente definido como 12 
nsel,s,loc,z,0 ! todos los nudos de la base 
d,all,uz,O ! con descenso impedido 
nsel,s,loc,z,19.70,19.80 ! todos los nudos de la 
cota de la linterna 

! con carga de linterna. 

! los extremos tienen la mitad 


f,a11,fz,-1.5 
nsel,r,loc,y,0,2,2 
f,all,fz,-0.75 


ANSYS 


nsel,all 

acel,0,0,1 Ipeso propio 

nrotat,1,315,1 ! preparo sistema de coordenadas 
locales en nudos paralelas a las 
cilindricas 

nsel,s,loc,y,0 ! 

d,all,uy,0 ! movimientos circunferenciales 

impedidos 
nsel,s,loc,y,2 ! 
d,all,uy,0 


¡RESISTENCIA FRENTE A EMPUJES 
¡Diferentes posibilidades 
11: mediante coacción infinita en la base, 
! dado un desplazamiento radial nulo 
12: mediante coacción en la base, dado un 
|. desplazamiento radial fijo no nulo 


! MODO 1: desplazamiento horizontal de apertura 
nulo 
' nsel,s,loc,z,0 
! d,alljux,0.00 
MODO 2: desplazamiento horizontal de apertura no 
nulo 
nsel,s,loc,z,0 
d,all,ux,0.01 
nsel,all 
/pbc,all,1 
eplot 


!lcm 


ISALVAR Y SALIR 
save 
finish ! salir del preprocesador 

¡PROCESADOR DE SOLUCION 

Imodulo de solucion estático, Newton Rapson 
/solution 
antype,static ! análisis estático 
outres,all,a]] 
neqit,40 
pred,on 
nropt,auto 
nsubst,15,20  ! especifica el numero y tamaño de los 
saltos de carga al estar activado los 
saltos automaticos (autots,on) 

! (15)tamaño del primer paso de carga, 
1 (20) numero de pasos máximo 
autots, on l activa saltos de carga automaticos 
solve 

finish 


MÉTODO DE LOS ELEMENTOS FINITOS 


44. CUPULA DE SAN PEDRO 11.1 


IPOSTPROCESADOR 
/postl 
set,l 
/pbc,al1,0 
esel,s,elem,,1,160,1  ! selecciona un subconjunto (s) 
de elementos por su numeración 
(elem), desde el 1 al 160 tomados 
delenl 
/view,1,0,1,0 
nsel,s,loc,z,-0.2,0.01 ! selecciona un subconjunto (s) 
de nodos de coordenada z entre — 
0.2 y 0.01 
treacciones en la base 
prrsol,fx  ! imprime las reacciones en los nudos en el 
sentido x, (fx) 
prrsol,fy  ! Idem, (fy) 
prrsol,fz  ! Idem, (fz) 
nsel,s,loc,y,0,0 
!tensiones. 
rsys,11 
pinsol,sx ! gráfica de resultados de las tensiones sx 
t en los nodos y elementos seleccionados 
! Idem, tensiones sy 
! Idem, tensiones sz 


pinsol,sy 

pinsol,sz 

nsel,all 

/device,vector,on ! opciones de dispositivo gráfico, 
l activa el modo VECTOR, las areas y 
volumenes se representan alambricas 

plerack ! representa las grietas 

/device,vector,off 


'ALGUNOS CALCULOS 
!tensiones en la base 


Ipath,1,22,43,64,85  ! define un camino de nodos 

para operaciones posteriores 

pdef,sig base-z,s,z  ! define e interpola entre los 
nodos del camino definido los valores de 
tensión en la dirección Z 

pdef,sig_cir,s,y ! Idem, tensión en la dirección Y 

plpath,sig_base-z,sig_cir 


ltrazado de tensiones sobre el eje de la cúpula 


Ipath,43,45,47,49,51,53,55,57,59,63 
pdef,s_z,s,z  ! tensión en la dirección Z 
pdef,s_y,s,y  ! tensión en la dirección Y 
plpath,s_z,S_y 
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'Trazado de tensiones en dirección Z en cinco cortes 
a lo largo de la cúpula. 
rsys,11 
dsys,11 
Ipath,1,3,5,7,9,11,13,15,17,21 
pdef,s_z1,s,z 
pdef,save  !evita que al definir un nuevo camino 
de nodos se pierdan los valores calculados 
Ipath,22,24,26,28,30,32,34,36,38,42 
pdef,s_z2,S,z 
pdef,save 
Ipath,43,45,47,49,51,53,55,57,59,63 
pdef,s_z3,5,z 
pdef,save 
Ipath,64,66,68,70,72,74,76,78,80,84 
pdef,s_z4,5,z 
pdef,save 
Ipath,85,87,89,91,93,95,97,99,101,105 
pdef,s_z5,s,z 
pdef,save 
pcalc,add,cero,s_z1,s_z1,1,-1 
pdef,save 


plpath,s_zl,s_z2,s_z3,s_z4,s_z5,cero 


! Mas interesante aun: trazado de espesores y 
! posicion del poligono de presiones: y grafica de 
tensiones N fi 


!Para espesor: distribución lineal de 3m. para s=0, y 
2m. para s=26.3 e=3-(s/26.3) 


pcalc,add,esp2,s,,-0.0190114,,1.5 
pcalc,add,jul,s_z2,5_z3,2,2 
pcalc,add,ju2,jul,s_z4,1,2 
pcalc,add,ju3,ju2,s_z1,1,1 
pcalc,add,ju4,ju3,s_z5,1,1 
pcalc,mult,ntotzju,ju4,esp2,0.25 
pcalc,add,jul,s_z1,s_z2,(-5/96),(-1/16) 
pcalc,add,ju2,s_z4,s_25,(1/16),(5/96) 
pcalc,add,ju3,jul,ju2,1,1 
pcalc,mult,ju4,esp2,esp2,4 
pcalc,mult,momju,ju4,ju3,1 
pcalc.div,excju,momju,ntotzju 
pcalc,add,-esp2,esp2,esp2,-.5,-.5 


plpath,esp2,cero,excju,-esp2 
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ESTRUCTURAS ll 


plpath,cero,ntotzju,momju 


'Trazado de tensiones en dirección Y en cinco cortes a 
lo largo de la cúpula. Y grafica de tensiones N teta 


rsys,11 
dsys,11 


ipath,1,3,5,7,9,11,13,15,17,21 

pdef,s_yl,s,y 

pdef,save l evita que al definir un nuevo camino 
de nodos se pierdan los valores calculados 
Ipath,22,24,26,28,30,32,34,36,38,42 

pdef,s_y2,S,y 

pdef,save 

Ipath,43,45,47,49,51,53,55,57,59,63 

pdef,s_y3,S,y 

pdef,save 

Ipath,64,66,68,70,72,74,76,78,80,84 

pdef,s_y4,s,y 

pdef,save l evita que al definir un nuevo camino 
de nodos se pierdan los valores calculados 
Ipath,85,87,89,91,93,95,97,99,101,105 
pdef,s_yS,s,y 

pdef,save 

pcalc,add,cero,s_y1,s_y1,1,-1 

pdef,save 

plpath,s_yl,s_y2,5_y3,s_y4,s_y5,cero 


pcalc,add,jul,s_y2,s_y3,2,2 
pcalc,add,ju2,jul,s_y4,1,2 
pcalc,add,ju3,ju2,s_y1,1,1 
pcalc,add,ju4,ju3,s_y5,1,1 
pcalc,mult,ntotyju,ju4,esp2,0.25 


plpath,cero,ntotyju 


!Trazado de tensiones en dirección XZ (tangenciales) 
en cinco cortes a lo largo de la cúpula. , grafica de 
tensiones tangenciales y suma de tangenciales con N 
fi, y todas juntas 


rsys,11 
dsys,11 
Ipath,1,3,5,7,9,11,13,15,17,21 


pdef,s_xz1,S,xz 
pdef,save ! evita que al definir un nuevo camino 


ANSYS 


de nodos se pierdan los valores calculados 
Ipath,22,24,26,28,30,32,34,36,38,42 
pdef,s_xz2,S,XZ 

pdef,save 

Ipath,43,45,47,49,51,53,55,57,59,63 
pdef,s_xz3,S,Xz 

pdef,save 

Ipath,64,66,68,70,72,74,76,78,80,84 
pdef,s_xz4,5,XZ 

pdef,save l evita que al definir un nuevo camino 


de nodos se pierdan los valores calculados 
Ipath,85,87,89,91,93,95,97,99,101,105 
pdef,s_xz5,S,xz 

pdef,save 

pcalc,add,cero,s_xz1,s_xz1,1,-1 

pdef,save 
plpath,s_xz1l,5_x22,5_xZ3,5_xZ4,5_xz5,cero 


pcalc,add,esp2,s,,-0.0190114,,1.5 
pcalc,add,jul,s_xz2,5_xz3,2,2 
pcalc,add,ju2,jul,s_xz4,1,2 
pcalc,add,ju3,ju2,s_xz1,1,1 
pcalc,add,ju4,ju3,s_xz5,1,1 
pcalc,mult,ntotxzju,ju4,esp2,0.25 
pcalc,add,ntotzxzju,ntotxzju,ntotzju,1,1 


plpath,cero,ntotxzju 
plpath,cero,ntotzxzju,ntotxzju,ntotzju,ntotyju 


finish 
FIN 


MÉTODO DE LOS ELEMENTOS FINITOS 


44. CUPULA DE SAN PEDRO 11.2 


MODO 1: 
SIN DESPLAZAMIENTOS HORIZONTALES EN LA BASE 
A | 


-309.598 
-266.295 
-222.992 
-179.69 
-136.387 
-93.084 
-49.782 
-6.479 
36.824 


A00ORCOSA 


GRAFICA DE TENSIONES SY (N teta) 


ANSY3 $S.3 
MAT 13 1999 
15:01:08 
POST: 
s7EP=1 

sua =S 
TIME=1 
PATH PLOT 
NCD:=35 
NOD2=1C05 


SAY PEDRO 


GRAFICAS DE NORMALES TOTALES Y MOMENTOS (N fi y mom) 


ESCUELA TÉCNICA SUPERIOR DE ARQUITECTURA 


ci DEPARTAMENTO DE ESTRUCTURAS DE EDIFICACIÓN 
U.P.M. MADRID 


ANSYS $.3 

MAY 13 1999 
15:56:13 

SODAL SOLUTION 


DMX =.001538 
SMN =-416.848 
3MX =-42,98 
-416.848 
-375.196 
-333.544 
-291.892 
-250.24 
-208.533 
-168.336 
-125.234 
-83.632 
-41.98 


EUDOSUOEN 


GRAFICA DE TENSIONES SZ (N fi) + 


ANSYS 5.3 
| may 12 1999 
15:59:55 
POST 

| STEP=1 

3U2 =S 
TIME-1 
PATH PLOT 
ROD1=85 

| sOD2=105 


3AN PEDRO 


GRAFICA DE EXCENTRICIDAD: POSICION DE LA RESULTANTE 


ESTRUCTURAS HI 


ANSTS 5.3 

MAY 13 1999 
19:57:13 

CRACKS AMD CRUSHIJG 
3T29=1 S 
308 =S 

TIME=1 


SAN PEDRO 


REPRESENTACION DE GRIETAS (ninguna en este caso) 


Í NOD2=105 


== 

¡ dIST=.7S 
xXx ..s 
Y -.s 
2 “Ss 
i 2-BUFF2R 


507.014 
674.463 


711.90 


SAN PEDRO ES 
GRAFICAS TOTALES DE: N fi, N teta, N tangencial y (N fi+ N tang) 


MÉTODO DE LOS ELEMENTOS FINITOS 
ANSYS 
A4. CUPULA DE SAN PEDRO 11.3 


MODO 2: 


CON DESPLAZAMIENTOS HORIZONTAL DE 1CM. EN LA BASE 


5.3 
MAT 13 1999 
15:37:08 
NODAL SOLUTION 
1 
| SUB =8 
TME=1 
(AVG) 
11 
DMZ =.010315 
=-1260 
=414.615 
-1280 
-1074 
-888.247 
702.124 
516.001 
¡ -329.878 
-143.755 
42.368 
228.492 
414.615 
SAN PEDRO 
GRAFICA DE TENSIONES SY (N teta) 
1s 5 
| MAY 23 1999 
15:35:55 
POST1 
STEP=1 
SUB =8 
ERAS TIME=1 
SN PATE PLOT 
6.0009 , KOD1=8S 
NOD2=105 
459.53: zY -=1 
DIST=.7S 
as. N, xr o -.S 
Ñ 6 TE o =.S 
5 MOeJu a 
sea 2-BUFTER 
Ñ 
3.072 a e 
Ñ F, MTOTZSG 
AN a E, j 
_ == 
2.97 ed 
587.98 E 
702.003 e 
/ 
236.511 
o 5.257 10.514 15. m.028 26.206 
2.028 13.143 14.400 23.037 
DIST 
SAN PEDRO 


GRAFICAS DE NORMALES TOTALES Y MOMENTOS (N fi y mom) 


U.P.M. 


DEPARTAMENTO DE ESTRUCTURAS DE EDIFICACIÓN 
ESCUELA TÉCNICA SUPERIOR DE ARQUITECTURA 


MADRID 


SAN PEDRO 


GRAFICA DE TENSIONES SZ (N fi) 


SAN PEDRO 


GRAFICA DE EXCENTRICIDAD: POSICION DE LA RESULTANTE 


ESTRUCTURAS Il 


ANSYS $.3 

MAT 13 1999 

15:39:23 

NODAL SOLCTION 
1 


32 AVG) 


DMX =.010315 

SMM =-926.185 

311.795 
-926.185 
-788.532 
-651.079 
-513.525 
-375.372 
-238.318 
-100.865 
36.688 
174.242 
311.795 


AODOBEOAN + 


AN3YS 5.2 


2Tr -. 
2-BUFFER 


SAN PEDRO 


NOD1=8S 


1025.10 
1196. 150 | 
1 
¡46.390 | NTOTYIU 
1 
 ————————— >= 
1536.43: 
o 5.2357 10.514 15.771 21.028 26.246 
2.028 7.889 13.143 18.400 1.097 
DIST 


E A A 
GRAFICAS TOTALES DE: N fi, N teta, N tangencial y (N fi+ N tang) 


MÉTODO DE LOS ELEMENTOS FINITOS 
ANSYS 
44. CUPULA DE SAN PEDRO 


Qrorar=107.602 KN SOLUCION 3:MODO 1 ANSYS 


SOLUCION 1: membrana a sin apertura en la base 


y | 197 
56.2 ( anillo linterna 
- A y anillos si los hay 


Qrorar=106.973 kN SOLUCION 3b: MODO 2 ANSYS 
SOLUCION 2: membrana apertura radial en la base 1cm 
DEPARTAMENTO DE ESTRUCTURAS DE EDIFICACIÓN : K 
LSE>  rscusLa récraca SUPERIOR DE ARQUITECTURA ESTRUCTURAS Ill o AS 
iS UPM MADRID ANSYS 


44. CUPULA DE SAN PEDRO 143 


¡PREPROCESADO, ENTRADA DE DATOS 
/prep7 
/window,1,full 
/view,1,0,0,1 
/vup,l,y 
/title, CANALIZACION DE AGUA 


!TIPO DE ELEMENTO 
et,1,4 ! BEAMA, admite grandes deformaciones 
| para analisis de 2$ orden 
r,1,300*300,300*300*300*300/12,81e8/12,300,300 
! AREA [mm2], lzz [mm4], Iyy [mm4], 
! Tkz [mm], Tky [mm] 
! toma como inercia a torsién la suma Izz+lyy 


IMATERIAL madera 
mp,ex,1,10 
mp,dens,0 ! no considera el peso propio 


¡NODOS 
n,1,0,0 
n,13,1000,12000 
n,16,4000,12000 
n,28,5000,0 
£111,1,13 ! genera una línea de nudos entre el 1 y el 13 
fi11,13,16 
f111,16,28 


¡ELEMENTOS 

type, 1 

real, 1 

84,2 

egen,12,1,1,1,1 
! generacion de elementos 
! 12 elementos generados (incluido el original), 
! incrementando la numeracién de los nodos en 1, 
i copiando del 1 al 1, tomados de 1 en uno 

e,13,14 

egen,3,1,13,13,1 

e,28,27 

egen,12,-1,16,16,1 

/pnum,node, 1 

/pnum,elem,1 

eplot 


¡CARGAS 
f£,13,£x,2 ! carga puntual en nodo 13, eje x, valor 2 [kN] 


f,13,fy,-160 
f,16,fx,2 
f,16,fy,-160 
/pbc,all, 1 


¡CONDICIONES DE APOYO 
nsel,s,loc,y,0 
d,alljux,0  ! articulaciones en los apoyos 
d,all,uy,0 
nsel,all 
d,all,uz,0 ! sin desplazamiento en Z, estructura arriostrada 


¡SALVAR Y SALIR 
save 
finish ! salir del preprocesador 


¡RESOLUCION 
! por la fuerte carga que soporta y su limitada rigidez, 
! debera analizarse a estabilidad. 
! Para estudiar la estabilidad hacemos un analisis 
! no lineal que llega al equilibrio cuando 
! en un proceso iterativo se converge. 
! La carga se introduce en 20 pasos. 
/solution 
antype,static 
nigeom,on l activa la no linealidad geometrica, 

! considera la posicion deformada de la estructura 
outres,all,all ! controla que la solución se escriba en la base de 
datos 

! (all) todos los casos de solución 

! (all) todos los saltos de carga 
neqit,50  !n” máximo de iteraciones de equilibrio= 50 
nsubst,20  !20 pasos de carga 
solve 
finish 
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12 CANALIZACIÓN DE AGUA 19.1 


CONJETURAS Diagramas de momentos 4.- MOMENTOS Y COMPROBACIÓN 
inicial sin giro con giro en los nudos 


1.- FUERZA DESESTABILIZADORA se La carga crítica de 2 kN producirá el diagrama: 
Ko 15 . en 1” orden en 2* orden 


$ 208-5 ki 
os acors 208.3 Kky 


Si la deformación es 0 


24.6 kS-m 


2.- FUERZA ESTABILIZADORA 
Deformamos el pórtico la misma Ó para ver la 
fuerza generada por su rigidez. 


El= 6750 kN.m* 
1 ? 
= —3=1405 a ] 
di 12* ? ' te 287 y las tensiones serán: 
18 27,4 + 1000 
M,=6El 7%=7508 AAA - id 
2 76 O, 03-300-300/6 6,07 N/mm” flexión 
h _ 160 » 1000 
3.- CARGA CRÍTICA Y AMPLIFICACIÓN á 300300 


Nes =28,58; N;= 342 kN (TOTAL 684) SroraL= 7,87 N/mn? 


= 1,8 N/mn? compresión 


factor amplificación: Coeficiente de seguridad de la estructura: 12/7,87 = 1,52 


A AA 
a= a 


342 
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¡SALIDA DE RESULTADOS 
/postl 

! Tablas de valores 
set,l ! último paso de carga 
etable,n1,smisc, 1 ! axiles en el nodo i 
etable,n2,smisc,7 ! axiles en el nodo j 
etable,t1,smisc,2 ! cortantes en el nodo i 
etable,t2,smisc,8 ! cortantes en el nodo j 
etable,m1,smisc,6 ' momentos en el nodo i 
etable,m2,smisc,12  ! momentos en el nodo j 
smult,ma,m!1,,-.001 ! cambio de unidades a [kN.m] 
smult,mb,m2,,-.001 

IGraficos 
pldisp,1 ! muestra la estructura deformada 

! y la original de referencia (1) 

*ask,p,¿sigo (gráfica de axiles)? 
plis.n1,n2  ! gráfica axiles [kN] 
*ask,p,¿sigo (gráfica de cortantes)? 
plis,t1,t2  ! gráfica cortantes [kN] 
*ask,p,¿sigo (gráfica de momentos)? 
plis,ma,mb  ! gráfica momentos [kN.m] 

Imprimir resultados 
prdisp ! imprime los desplazamientos 
pretab,n1,n2,t1,t2,ma,mb 
! imprime los valores de las tablas relacionadas 
! imprime las reacciones en los nudos (f): fuerzas 
! imprime las reacciones por elementos 

! (£): suma de las fuerzas nodales 

*ask,p,¿sigo (evolución de gráfica de momentos)? 
*do,ciclo,1,20 
I*ask,p,¿sigo? 
set,1,ciclo 
etable,m1,smisc,6 
etable,m2,smisc,12 
smult,ma,m1,,-.001 
smult,mb,m2,,-.001 
/contour,1,9,-40,,40 
*cet,escalal,elem,12,smisc, 12 
*set,escl,((2)*escala1/(40000)) 
plis,ma,mb,esc] 
*enddo 
'finish 


prrsol,f 
presol,f 


! cambio de unidades a [kN.m] 


/post26 ! entra en el posp. 26 (time-history) 
nsol,2,13,u,x,desp-hor !define var 2 como ux nodo 13 


plvar,2 !gráfica ux nodo 13 en el proceso de carga 
finish 
rs... DEPARTAMENTO DE ESTRUCTURAS DE EDIFICACIÓN 
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ESFUERZOS Y DEFORMACIONES 


Primer orden 


HI — +35 25 57 —— - 249% 
Fu AO MUA 23. SY 
¿ : f 
by a 
EQUILIBRIO EN NUDOS 
Primer orden Segundo orden 
a AS 96 KS dE Ape KS 
a q ECN EA ARA na SS 
q . 
An Y $ 15,5 4, + al? 44? 
3 150,4 (69,6 ! ds las ( Ñ 
(AO part 
104 rs ua W 1355 o e A + 
Leda M Guslunidas 4 
Mo X Y ES 4 20 
Log 18 Luo as 0 Lor 


ESTRUCTURAS III 


Segundo orden 


MÉTODO DE LOS ELEMENTOS FINITOS 
ANSYS 
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12.3 


'PREPROCESADO, ENTRADA DE DATOS 
! búsqueda iterativa de la forma de una lámina 
/prep7 
/window, 1, full 
/view,1,2,-3,5 
/vup.!1.z 
/title, PLACA CUADRADA 


!TIPO DE ELEMENTO 
et,1.63 ! elemento tipo shell63, plano de 4 nodos 
r.1,7 ! espesor [cm] 


¡MATERIAL hormigón 
mp.ex.1.2.1e3 ! modulo de elasticidad [kN/cm2] 
mp.nuxy,1..] ! coef. de poisson 


mp.dens.1,2.5e-5  ! densidad [kN/cm3] 


¡DEFINICIONES 
an=11*11 ' numero de nudos 
nciclos=20 ! numero de ciclos 

! matrices para las coordenadas de los nudos 
*dim,x,array,nn t define la matriz x de nn 
*dim,y.array,nn ! define la matriz y de nn 
*dim,z,array,nn ! define la matriz z de nn 

! matrices para los desplazamientos de los nudos 
*dim,dl,array,nciclos+]1 
*dim.d2.array,nciclos+1 
*dim.d3,array.nciclos+1 


¡GEOMETRIA 

n.1,0,0.0 ! define el nodo 1 de coordenadas (0,0,0)cm 
n.11,.3000,0,0  ! define el nodo 10 de coordenadas (3000.0.0)cm 
£1,1.11 ! genera una línea de nodos entre 1 y 10 
ngen.11,11.1,11.1,0,300,0,0 

! ngen: genera nodos según un patrón: 

! 11 conjuntos de nodos(incluido el original), 

! incrementando la numeración en 11, 

! copiando del nodo 1 al 11 tomados cada 1, 

' incrementando X en 0, Y en 333.3cm y Z en 0 
nplot ! dibuja los nodos en pantalla 


¡GUARDA COORDENADAS INICIALES 

*do.n,1,nn, 1 ! bucle de n=1 hasta nn, incremento de 1 
x(n)=nx(n) ! las funciones nx(n). ny(n), nz(n), 
y(n)=ny(n) ! devuelven las coordenadas x. y. z del nodo n 
z(n)=nz(n) 

*enddo 


DEPARTAMENTO DE ESTRUCTURAS DE EDIFICACIÓN 


¡GENERACION DE ELEMENTOS 
type.! 
real. 1 
€e,1,2,13 
e,1,13,12 
egen.5.1.1,2,1 
! generación 5 conjuntos de elementos 
! (incluido el original), incrementando la 
! numeración de los nodos en 1, 
! copiando del 1 al 2, tomados de 1 en uno 
egen.5,11,1,10,1 
egen.2,60.1.50,1 
6:6.7,17 
e,7.18,17 
egen.5.1,101,102,1 
egen.5,11,101,110,1 
egen.2,50,101.150,1 
/pnum,node,0  ! desactiva la numeración de los nodos 
/pnum.elem.!  ! activa la numeración de elementos 
eplot ! dibuja elementos 


! CONDICIONES DE APOYO 
/pbc.u.l - L activa los símbolos de coacciones 
! a desplazamiento 
nsel.s.node,. 1 
nsel.a.node.,.1 1 
nsel.a.node.,111 
nsel.a,node..121 
d,all.uz.0 
d,all,ux.0 
d.all,uy.0 
nsel.all 


¡CARGAS 
acel.0,0.-1 ! asigna una aceleración -1 en el sentido del eje Z 
! se pone la carga hacia arriba para ver la lamina 
! invertida 
eplot 


ISALVAR Y SALIR 
save 
finish ! salir del preprocesador 


¡BUSQUEDA DE LA FORMA: PROCESO ITERATIVO 
/input.busqueda.dat 


¡NCOOR.DAT 
¡ASIGNA LAS COORDENADAS DE 1.0S NUDOS 
*do.i.1.nn.1l ! bucle de i=1 hasta numero de nodos. 
l incrementándolo en 1 
n.1.x(1).y(1).z(1) ! define las coordenadas de los nodos 
*enddo ! fin del bucle 


! CAMBIO.DAT: CAMBIO DE COORDENADAS 
! cambia los valores de las matrices x, y. z. sumando al 
! valor existente el desplazamiento * factor de escala. 
llas funciones ux(n), uy(n). uz(n) 
! devuelven los desplazamientos del nodo n 
*do.n.1,nn.l ! bucle de n=1 hasta nn. incremento de 1 
x(n)=x(n)+fe*ux(n) 
y(n)=y(n)+fe*uy(n) 
z(n)=z(n)+fe*uz(n) 
*enddo 


! BUSQUEDA.DAT PROCESO DE BUSQUEDA DE LA FORMA 
! 1? ITERACION SOLUCION Y CAMBIO DE COORDENADAS 
/solution 
solve ! ordena resolver el sistema 
finish ! salir del procesador 
fe=.2 ! factor de escala 
/input.cambio.dat  ! cambia los valores de las matrices x. y. z 
/prep7 
/input.ncoor.dat ! lee el fichero ncoor.dat 

! que asigna las coordenadas de los nudos 
/pnum.node.0 l activa la numeración de nodos 
eplot ! dibuja los nodos en pantalla 
finish 


![TERACIONES SIGUIENTES: 
Incilcos SOLUCIONES Y CAMBIOS 
fe=.5 ! factor de escala 
*do.nc.1.nciclos. | 
/solution 
solve l ordena resolver el sistema 
finish ! salir del procesador 
dlí(nc)=uz(61) ! guarda la deformación del nodo 61 
d2(nc)=uz(6) ! guarda la deformación del nodo 6 
d3(nc)=uz(3) ! guarda la deformación del nodo 3 
/input.cambio.dat ! cambia las matrices x. y. z 
/prep7 
/input.ncoor.dat  ! lee el fichero ncoor.dat 
eplot ! dibuja los elementos en pantalla 
finish 


LA ESCUELA TÉCNICA SUPERIOR DE ARQUITECTURA ESTRUCTURAS TI METODO DE LOS ELEMENTOS FINITOS 
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13. BÚSQUEDA DE FORMAS EN LAMINAS 13: 1 


¡"ANALISIS DE LA ESTRUCTURA RESULTANTE ! limprimir 
/solution pretab.s1.s2  ! imprime los datos de las tablas TENSIONES PRINCIPALES: S1 TKN/omA 
solve ! ordena resolver el sistema I sl. s2 [kN/cm2] 
finish ! salir del procesador prrsol.f | imprime las reacciones en los nudos. (f) [kN] 
dl(nciclos+1)=uz(61) ! guarda la deformación del nodo 61 prnsol.u ! imprime desplazamientos (u) en los nodos [cm] 
d2(nciclos+1)=uz(6) | guarda la deformación del nodo 6 *status.d | ¡ imprime los valores de la matriz dl donde se 
d3(nciclos+1)=uz(3) ! guarda la deformación del nodo 3 ¡ guardan las deformaciones del nodo 61 
*status.d2 ¡ ídem nodo 6 
!'POSTPROCESADOR *status.d3 ¡ ídem nodo 3 
I¡SALIDA E INTERPRETACION DE RESULTADOS finish 
¡ESTRUCTURA 1 
/postl 
shell.bot l especifica la localización en que se dan los 
! resultados de tensiones en el elemento laminar 
! opciones: (top) cara superior.(mid) medio x $e 
! (bot) cara inferior S Es ES L3SE TA 
etable.sxx.s,x  ! prepara la tabla denominada sxx . y 200Sa 
! con las tensiones (s). componente X 
l estas tensiones están en coordenadas globales 
etable.sy y.S.y 
etable.Sxy.S.Xy 
etable.sl.s.l ! prepara la tabla denominada sl 


-.133E-11 
.031883 
063767 
09565 

. 127534 

. 159417 
«191301 
223184 
.258068 


ONIAAAnaD * 


' : 
' con las tensiones (s), componente 1 EN 
í tensión principal 


etable.s2.s.2 
set.1.1 TENSIONES PRINCIPALES: S2 LN lem 


lgraficas 
/pbc.u. 1 ! pone los símbolos de las coacciones a 
' desplazamientos (u), en activo(1) 
/dscale. 1 ¡ escala de la representación de la deformada 
¡ en pantalla DEFORMACION EN LOS NODOS SELECCIONADOS 
pldisp.2 ! muestra la estructura deformada ceca ps e 
! y el contomo de la original (2) 
pletab.sxx.avg ! gráfica de los datos de la tabla sxx, 
! (avg): tomando en el interior del elemento un 
! valor como promedio lineal 
pletab,syy.avg 
pletab.sxy,avg 
pletab.sl.avg  ! gráfica de los datos de la tabla sl. 
! (avg): tomando en el interior del elemento un 
' valor como promedio lineal 
pletab.s2.avg 
/vscale. 1.3 l escala la representación de vectores en pantalla 
plvect.s...,..on — ! gráfica de las tensiones principales como vector 


deformación 


-.006117 
-.00255 
. 001016 
004582 
.008148 
«011715 
.015281 
.018847 
022414 


Iteraciones 


OOONOA Ma, * 
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u) 
(59) 


¡PREPROCESADO, ENTRADA DE DATOS 
TODOS LOS ELEMENTOS SON CUADRANGULARES 


/prep7 l entra en el preprocesador 
/window.1.full 

/view.1.2.-3.5 

/vup.T.z 


hítle, PLACA CUADRADA 


TIPO DE ELEMENTO 


et.1,63 l elemento tipo shell63, 
! plano de 4 nodos 
mui ! espesor [cm] 
¡MATERIAL hormigón 
mp.ex,1.2.1e3 ! modulo de elasticidad [kN/cm2] 


mp.nuxy.!..l 
mp.dens.1.2.5e-5 


! coef. de poisson 
! densidad [kN/cm3] 


¡GEOMETRIA 

n,1.0.0,0 ! define el nodo 1 de coordenadas (0,0.0)cm 
n.11.3000,0.0  ! define el nodo 10 de coordenadas (3000.0.0)cm 
£11,1,.11 ! genera una línea de nodos entre 1 y 10 
ngen.11.11,1.11.1,0,300,0,0 

! ngen: genera nodos seg£n un patron: 

! 1! conjuntos de nodos(incluido el original). 

! incrementando la numeración en 11. 

! copiando del nodo | al 11 tomados cada 1. 

! incrementando X en (0), Y en 333. 3cmyZen0 


nplot ! dibuja los nodos en pantalla 
"DEFINICIONES 
nn=11*11 


! define la matriz x de nn 
! define la matriz y de nn 
! define la matriz z de nn 


*dim,x.array.nn 
*dim,y.array.nn 
*dim.z.array,nn 


'GUARDA COORDENADAS 
*do.n.1.nn.! ! bucle de n=1 hasta nn. incremento de 1 
x(n)=nx(n) 
y(n)=ny(n) 
z(n)=nz(n) 
*enddo 


generacion de elementos 


type. 1 
real,l 
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e.1,2,13 
e.1,13,12 
egen.5,1.1.2,1 
! generación 5 conjuntos de elementos 
! (incluido el original). incrementando la 
! numeración de los nodos en 1. 
! copiandce, del | al 2. tomados de | en uno 
egen.5.11.1,10,1 
egen.2.60,1,50,1 
e.6,7,17 
e.7,18,17 
egen.5.1,101,102,1 
egen.5.11.101.110,1 
egen.2,50,101,150,1 
/pnum.node.0 ! desactiva la numeración de los nodos 
/pnum.elem. |! l activa la numeración de elementos 
eplot ! dibuja elementos 


! CONDICIONES DE APOYO 
/pbc.u.] ! muestra los símbolos de coacciones 
¡ a desplazamiento 
nsel.s.node..! 
nsel.a.node..11 
nsel.a.node..111 
nsel,a.node..121 


d.all.uz.0 

d.all,ux.0 

d.all,uy.0 

nsel.ali 

¡CARGAS 

acel,0,0.-1 ! asigna una aceleración -.!| en el sentido del eje Z 
! se pone la carga hacia arriba para ver la lamina 
! invertida 

eplot 


ISALVAR Y SALIR 
save 
finish ! salir del preprocesador 


' GUARDA LAS COORDENADAS SOLUCION 0: INICIAL 
*dim.x0,array.nn  ! define la matriz x0 de nn 
*dim.y0.array.nn  ! define la matriz yO de nn 
*dim.z0,array.nn  ! define la matriz z0 de nn 


*do.n.1.nn, | ! bucle de n=1 hasta nn. incremento de 1 


ESTRUCTURAS III 


x0(n)=x(n) 
y0(n)=y(n) 
z0(n)=z(n) 
*enddo 
¡ANALISIS 1 
/solution ! entra en el procesador de solución 
solve ! ordena resolver el sistema 
finish ! salir del procesador 


!VER DEFORMADA 
/postl ! entra en el postprocesador 


/dscale.1..2 l escala la representación de la deformada 
pldisp.2 ! muestra la estructura deformada 
finish 


¡GUARDA DEFORMADA 1 
*dim.xl.array.nn  ! define la matriz xl de nn 
*dim.yl.array.nn  ! define la matriz y! de nn 
*dim.zl.array.nn  ! define la matriz z1 de nn 
fe=.2 ! factor de escala 
*do.n.!.nn.! ! bucle de n=1 hasta nn. incremento de | 
x1(n)=x(n)+fe*ux(n) 
y Hn)=y(n)+fe*uy(n) 
z1(n)=z(n)+fe*uz(n) 
*enddo 


!MATERIAL hormigon ligero 


/prep7 ! entra en el preprocesador 
/pnum.node. | l activa la numeración de los nodos 
nplot ! dibuja los nodos 

mp.ex.2,2.1e3 ! modulo de elasticidad [kN/cm2] 
mp.nuxy.2..l ! coef. de poisson 

mp.dens.2.0.0 ! densidad [kN/cm3] sin peso propio 


! ELEMENTOS NUEVOS: BARRAS 
et.2,4 l elemento tipo 2 BEAMA4 
1:2.900:2.787.2.107.1,1 
! constantes para la sección 2: AREA [cm2] 
!Izz [cm4]. lyy [cm4]. Tkz [cm]. Tky [em] 


! toma como inercia a torsión la suma lz7+ly 


Igeneracion de elementos 


mat.2 ! asigna el material 2 a los siguientes elementos 
type.2 ! asigna el elemento tipo 2 a los siguientes 

¡ elementos 
real.2 ! asigna las constantes 2 a los siguientes elementos 
e.1.13 ! define el elemento | que va del nodo ! al 17 


MÉTODO DE LOS ELEMENTOS FINITOS 
ANSYS 
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y 


14, 


egen.10.12.201.201.1 
! generación de 10 conjuntos de elementos. 
! incrementando los nodos en 12, 
' tomando del elemnto 201 al 201 cada | 

e.11.21 

egen.10,10,211.211.1 

eplot 

¡SALVAR Y SALIR 
save 
finish i salir del preprocesador 


¡ANALISIS 2 


/solution 
solve ! ordena resolver el sistema 
finish ! salir del procesador 


'VER DEFORMADA 


/postl 

/dscale, 1.4 ' escala la representación de la deformada 
pldisp.2 ! muestra la estructura deformada 

finish 


¡GUARDA DEFORMADA 2 
*dim.x2.array.nn  ! define la matriz x1 de nn 
*dim.y2.array.nn — ! define la matriz yl de nn 
*dim.z2.array.nn  ! define la matriz z1 de nn 
fe=1 ! factor de escala 
*do.n.1.nn, 1 t bucle de n=1 hasta nn. incremento de 1 
x2(n)=x(n)+fe*ux(n) 
y2(n)=y(n)+fe*uy(n) 
z2(n)=z(n)+fe*uz(n) 
*enddo 


! ELIMINA ELEMENTOS BARRA 
/prep7 
eplot 
edelc.201.220.1 


! COMBINA LAS COORDENADAS SOLUCION 1 Y 2 


fel=.5 ! factor de escala caso 1 

fe2=2 ! factor de escala caso 2 

*do.n.1,nn.] ! bucle de n=1 hasta nn. incremento de 1 
z(n)=z1(n)*fe1+z2(n)*fe2 

*enddo 
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¡ASIGNA LAS COORDENADAS DE LOS NUDOS 


*do.¡.1.nn.l ! bucle de i=1 hasta numero de nodos. 
! incrementandolo en | 
n..x(1).y(1).z(1) ' define las coordenadas de los nados 
*enddo ! fin del bucle 
eplot 
finish ! salir de preprocesador 


¡ANALISIS DE LA ESTRUCTURA RESULTANTE 3 


/solution 

solve ordena resolver el sistema 

finish ! salir del procesador 
¡POSTPROCESADOR 


ISALIDA E INTERPRETACION DE RESULTADOS 
¡ESTRUCTURA 1 

/post! 

shell.top l especifica la localizacion en que se dan los 


! resultados de tensiones en el elemento laminar 


! opciones: (top) cara superior.(mid) medio 
! (bot) cara inferior 
etable.s!.s.] ! prepara la tabla denominada s! 
! con las tensiones (s). componente 1 
! tension principal 


etable.s2.s.2 
set. 1.1 
Igraficas 
/pbc.u. | ! pone los simbolos de las coacciones a 
! desplazamientos (u). en activo(1) 
/dscale. 1 
pldisp.2 | muestra la estructura deformada 


! y el contorno de la original (2) 

pletab.sl.avg ! gráfica de los datos de la tabla s!. 
! (avg): tomando en el interior del elemento un 
! valor como promedio lineal 

pletab.s2.avg 


limprimir 
pretab.sl.s2 ! imprime los datos de las tablas 
sI, s2 [kN/cm2] 
prdisp ! imprime los valores de la deformacion 
prrsol.f ! imprime las reacciones en los nudos. (f) [kN] 
finish 


ESTRUCTURAS III 


MÉTODO DE LOS ELEMENTOS FINITOS 
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